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Abstract. This paper is on magnetic Schrodinger operators in two dimensional domains with 
corners. Semiclassical formulas are obtained for the sum and number of eigenvalues. The 
obtained results extend former formulas for smooth domains in 110) to piecewise smooth 
domains. 

1. Introduction 

The spectral analysis of magnetic Schrodinger operators in domains with boundary has been 
the subject of many research papers in the last two decades. Apart from the mathematical 
interest behind the study of their spectra, magnetic Schrodinger operators in interior /exterior 
domains with various boundary conditions arise in several models of condensed matter physics, 
as superconductivity |14 | l9l 116] . liquid crystals \12\ I18| and Fermi gases [6]. 

The present paper is devoted to the study of magnetic Schrodinger operators in domains with 
corners (piecewise smooth domains). The presence of corners in the domain has a strong effect 
on the spectrum of the operator. In particular, it is shown in [2j [13] that the presence of corners 
decreases the value of the ground state energy of the operator compared with the case of smooth 
domains. Discussion of this effect in the framework of superconductivity is given in [3|. 

We give in a simple particular case, a brief presentation of the semiclassical results proved 
in this paper. Suppose for simplicity that ft is a simply connected bounded domain in M 2 and 
that A is a vector field such that b = curl A is constant. Let Ph,n = — (hV — iA) 2 be the 
magnetic Laplacian in L 2 (Q) with magnetic Neumann boundary condition, N(Xh) and E(Xh) 
be the number and sum of negative eigenvalues of Ph,n — Xld. If the boundary of is smooth, 
it is proved in |10| [TT] that, as h — > 0, 

N(Xh) ~ci(6,A)|50|/i~ 1/2 , E(Xh) ~ c 2 (b,X)\dn\h 1 / 2 , (1.1) 

where c\ and C2 are two explicit constants. The formula for the number N(Xh) is valid for all 
A < 6, while that for the energy E{Xh) is valid for all X < b. 

In Theorem 14.11 of this paper, we prove that (jl.ip is still holding true when the domain Q 
is only piecewise smooth. The key to prove this extension is a rough estimate on the number 
of eigenvalues of a Schrodinger operator with constant magnetic field in a sectorial domain (see 
Lemma 12. 5p . 

The result of Theorem 14.11 shows that corners only affect the low-lying eigenvalues of the 
operator. As we will see in Theorem 13.11 corners create a few additional eigenvalues compared 
with smooth domains. 

As a consequence, we may say that in the semi-classical limit and the regime considered in this 
paper, the number of eigenvalues of magnetic Schrodinger operators for smooth and non-smooth 
domains are asymptotically the same. 
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The paper is organized in the following way. In Section [2j we collect some key results that 
will be used throughout the paper. Section [3] is about the number of low- lying eigenvalues of the 
operator in a domain with corners (here the regularity assumption on the domain is precisely 
stated). Finally Section 2] contains the semi-classical analysis that extends (jl.ip to domains with 
corners. 

Throughout the paper, the notation N(X, P) will be often used to denote the number of 
eigenvalues (counting multiplicities) of the operator P that are below A. 

2. Preliminaries 

2.1. Variational principles. Let H be a self-adjoint operator in a Hilbert space T~L (of domain 
D{H)) such that 

inf a css (H) > 
i/l(_ OOi 0)(-ff) is trace class . 
We shall need the following two simple variational principles concerning the operator H, which 
are frequently used in \17\ [8] . 



(H) 



Lemma 2.1. Let 7 be a bounded operator such that < 7 < 1 (in the sense of quadratic forms) 
and the operator H 7 is trace class. Then it holds that, 

tr^l^o)^)) <tr(ff 7 ). 

Lemma 2.2. Assume that the operator H satisfies the hypothesis (H). Then it holds that, 

N 

tr (#1(^0) (ff)) = inf J2(fj > H fj) > 

where the infimum is taken over all orthonormal families {fx, /2, . . . , /at} C D(H) and N > 1. 

2.2. A family of one-dimensional differential operators. Let us recall the main results ob- 
tained in [3[T3] concerning the family of harmonic oscillators with Neumann boundary condition. 
Given ( 6 8, we define the quadratic form, 

B 1 (R + )3u^q[Z](u)= [ \u'(t)\ 2 + \(t-£)u(t)\ 2 dt, (2.1) 

Jr + 

where, for a positive integer fc£N and a given interval ICI, the space B k {I) is defined by : 

B k (I) = {u G H k (I); t j u(t) G L 2 (I), Vj = 1, . . . , k}. (2.2) 

Since the quadratic form f)2. 1|) is closed and symmetric it defines a unique self-adjoint operator 
£[£]. This operator has domain, 

D(£[£]) = {ue B 2 (R + ); u'(0) = 0}, 

and is the realization of the differential operator, 

m = -d 2 + (t-o\ (2.3) 

on the given domain. We denote by {/ij(£)}j=a the increasing sequence of eigenvalues of £[£], 
which are all simple. By the min-max principle, we have, 

Mi(0= -f (2 . 4) 
ueB^+Wo \\u\\ L 2 {R+) 

It follows from analytic perturbation theory (see [7]) that the functions 
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are analytic. Furthermore, //i(0) = 1 and — 1| decays like exp(— £ 2 ) as £ — > +00 (see [1]) 
thus yielding that 

/ - 1) d£ = - / - 1]_ is finite. 

JO JR 

We define the constant : 

e = inW£). (2.5) 



In [7J, it is proved that 0o = /Ui(£o)> that £0 is the unique value at which the minimum Go is 
attained and that /u"(£o) > 0. 

An important consequence of standard Sturm-Liouville theory is recalled below (c.f. 
Lemma 2.1]). 

Lemma 2.3. The second eigenvalue of £[£] satisfies, 

inf fi 2 (C) > I- 



Notice that part of this conclusion is a consequence of the analysis of Dauge-Helffer [7J, who 
show that the infimum of [i2(Q is attained for a unique £2 & ^- 

2.3. Rough energy bound for the cylinder. Let us consider the operator 

Pn = -(W-^A ) 2 mL 2 (n), 

with 

Q = [0,5] x (0,/i 1/2 T). 

Functions in the domain of Ph,n satisfy Neumann condition at t = 0, periodic conditions at 
s £ {0, S} and Dirichlet condition at t = h l l 2 T. We assume that the vector field Ao is given by 

A (xi,x 2 ) = (-x 2 ,0) ■ 

In this particular case, the operator has compact resolvent, hence the spectrum consists of an 
increasing sequence of eigenvalues (ej)j>\ converging to +00. Note that the terms of the sequence 
(e,) are listed with multiplicities counted. Given A £ K, the energy 

£(X,b,S,T) =J2[hb(l + X) - (2.6) 

3 

is finite. The number of eigenvalues below (1 + X)hb is also of particular interest : 

jV{X,b,S,T) = Card{j : ej < hb(l + A)} . (2.7) 
In the lemma below, we state estimates of the energy £(X, b, S, T) and the number jV(X, b, S, T) . 

Lemma 2.4. There exist positive constants Tq and Xq such that, for all S > 0, b > 0, T > VBTq 
and X G (0, Ao], we have, 

£(X, b, S, T) < (1 + X)hb + l) , (2.8) 

ST 

M(X,b,S,T)<— - + 1. (2.9) 

The lemma is proved in |10| . although the estimate (|2.9p is not stated explicitly. Actually, 
(|2.8p is proved as follows. First (|2.9p is established by separation of variables and the variational 
min-max principle then the energy £ (A, b, S, T) is easily estimated as j\f(X, b, S, T) x (1 + X)bh, 
from which (12.81) follows. 
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2.4. Rough bounds for the sector. Let a G (0, 2tt), R> 0, h> 0, b> and 

ttR,h,a = {(rcos0,rsin0) G M 2 : < 9 < a, 0<r<h l/2 R}. (2.10) 

Consider the self-adjoint operator P^a = — (/iV — ibA^) 2 in L 2 (J7R/j a ). Functions in the 
domain of Ph, a satisfy Dirichlet condition on the boundary r = R and Neumann condition 

v ■ (/iV - ibA Q )u = 

on the boundary defined by 9 = or 9 = a. Here u is the unit outward normal vector on the 
boundary dQ,R,h,a an d Ao is the magnetic potential introduced in (|3.2p . 

The operator Ph, a has compact resolvent and its spectrum is discrete and consists of isolated 
eigenvalues (e_,) counted with multiplicities. Let A £ M and define, 

£corn(a, b, R,\) = ^2 [hb(l + A) - ej] + , (2.11) 
j 

and 

M corn (a,b,R,X) = Card{j : ej < hb(l + A) } . (2.12) 
We give rough estimates of £ COTn (a, b, R, A) and JV C oTn(ct, b, R, A) in the next lemma. 

Lemma 2.5. Given b > and a G [0,27r) ; there exist positive constants C, ho, Rq and X\ such 
that, for all R > Rg, h G (0, ho), and A G (— oo, Ai], we have, 

M corri {a, b, R,X)<C (R 2 + 1) , (2.13) 

and 

£ COII1 (a, b, R, A) < Chb (R 2 + l) . (2.14) 

Proof. Observe that the upper bound (|2.14p on the energy £" corn (a, b, R, A) follows immediately 
from the definition of £ corn (a, b, R, A) and the upper bound in (|2 . 13)) . Therefore, we establish the 
bound in (|2.13p . The method that will be used is introduced in [5] and based on a decomposition 
of the operator via a partition of unity and on the min-max variational principle. 

To simplify the notation, we denote by O and P the set £lR t h,a introduced in (|2.10p and the 
operator Ph, a respectively. 

Let Aq be a positive real number whose choice will be specified later. Suppose R > Aq. We 
cover O by two sets, 

A h l l 2 

U x (h) = {x G ft : |z| < A /i 1/2 }, U 2 (h) = {x G ft : — < \x\ < Rh 1/2 } . 

Consider a partition of unity of Q, 

2 2 c 

Y J x){x) = l and ^|V X ,(x)| 2 <-^inf], 

j=l j=l o 

where C\ > is a universal constant and suppxj C Ujfo 1 / 2 ). 

Let u G H 1 ^). The following decomposition formula holds true, 

2 

[ \{V - ibA )u\ 2 dx = V f (\(V -ibA ) Xj u\ 2 -h 2 \V X j\ 2 \u\ 2 )dx. 
Jn - =1 Jn 

This decomposition yields the inequality stated below by using the upper bound on |Vx?| : 

I \(V - ibA )u\ 2 dx > V f (\(V -ibA ) Xj u\ 2 -CiAQ 2 h\u\ 2 )dx. 
Jn . j Jn 
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By using the method in [5] and the variational min-max principle, we get that, 

Kom(a,b,R,X) < N(Ah,P Ul(h) ) + N(Ah,P U2(h) ), (2.15) 

where A = (1 + A)6 + Ci^q 2 , P^fh) = ~(hV -ibA ) 2 is the operator in L 2 (Ui(h)) with Dirichlet 
condition on r = /i 1 / 2 and Neumann condition on the other parts of the boundary of U\(h), and 
Pjj 2 {h) = —(hV — ibAo) 2 is the operator in L 2 {U 2 {h)) with Neumann boundary condition on 
9 = and 9 = a and Dirichlet condition on the other parts of the boundary. Notice that we use 
polar coordinates (r, 9) in the definition of the domains. 

Let Ao be as in the statement of Lemma 12.41 and select Ai G (0,Ao). In light of (|2.15p . 
Theorem 12.51 will be proved once the statement that follows below is shown to be true: 

3 C > 0, N(Ah,P Ul{h) ) < C and (2.16) 
N(Ah,P U2(h) )<C{R 2 + l) . (2.17) 

Proof of (|2.16p . Applying the re-scaling y = h~ 1 / 2 x changes the domain U\{h) to U\{Aq) and 
the operator Pjj^h) to the operator in L 2 (Ui(Aq)), 

Pu^Ao) = -(V-ibA ) 2 . 

More precisely, a number e is in the spectrum of Pu^h) if an d only if e/h is in the spectrum of 
PjjxlAo)- Thus, the term we wish to estimate is given as follows, 

N(Ah,P Ul{h) ) = N{A,P Ul{Ao) ). 

Since the operator Pjj^i) has compact resolvent and the number A is bounded independently of 
h G (0, 1], it follows that A^(A, P^m) is bounded as h varies in (0, 1] too. This establishes the 
upper bound in (|2.16p . 

Proof of (|2.17|) . As we shall see, the upper bound in (|2.17p follows mainly from Lemma [2.41 Let 
e be a positive real number such that, 



/an a 
0<e<wm[-,--- 



Recall the definition of the set 



u 2 (h) = {xen : < M < Rhl/2 } ■ 



We cover U2(h) by four sets, 



Vtih) = {(r cos 9,r sin0) G U 2 (h) : < 6 < -} 



a a 



V 2 (h) = {(r cos 9,r sin0) G U 2 (h) : j-e<6<- + e} 

a . 3a 



V 3 (h) = {(r cos 9, r sin 9) eU 2 (h) : - -e<6<— + e} 

3a 



Vi(h) ={(rcos0,rsin0) G U 2 {h) : — < 9 < a} . 

Define four operators P Vl (h), Pv 2 (h), p V 3 {h) and P Vi {h) m ^ 2 (^i(^)) 5 L 2 (V 2 (h)), L 2 (V 3 (h)) and 
L 2 (V4(/i)) respectively. The four operators are self-adjoint realizations of the differential operator 
— (/iV — ibAo) 2 . Functions in the domain of Py 1 (h) satisfy Neumann condition on 9 = and 
Dirichlet condition on the other parts of the boundary. Similarly, functions in the domain of 
Pv 4 (h) satisfy Neumann condition on 9 = a and Dirichlet condition on the other parts of the 
boundary. Functions in the domains of Pv 2 (h) an d Pv 3 (h) satisfy Dirichelt boundary condition. 
Thus, the bottoms of the spectra of the operators Pv 2 (h) an d Pv^ih) are larger or equal to bh. 
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Notice that the operators Py 1 (h) an d Pvi(h) are unitary equivalent and hence have same spectra. 
Also, the operators Pv 2 (h) an d -fV 3 (h) are unitary equivalent and have same spectra. 

We apply an argument similar to the one we did to obtain (|2.15p . By introducing a partition 
of unity supported in V±(h), V 2 (h), V^{h) and V^h), using the IMS decomposition formula and 
the variational min-max principle, we get a constant C 2 > that is allowed to depend on e but 
not on h and such that, 

4 

N(Ah,P U2{h) ) <Y,N(Ah,P Vj(h) ) 
2 

= 2j2N(Ah,Fv. (h) ), (2.18) 
i=i 

with A = A + CiAq 2 + C 2 h. 

Recall the number M(X,b, S,T) introduced in (|2.7p . This number counts the eigenvalues of 
the operator P v in the cylinder, where U = [0, S] x (0,h^ 2 T). 

Since the angle a/4 is acute and Dirichlet condition is imposed in the form domain of Py 1 ^ , 
functions in the form domain of Py 1 (h) can be extended by zero to functions defined in the 
cylinder, 

U = [0,Rh^ 2 ] x (0,i?/i 1/2 ). 

The extended functions are in the form domain of the operator Pjj. Thus, by the variational 
min-max principle, it is easy to see that, 

N(Ah,P Vl(h) ) < AA(A,6,i?/i 1 / 2 , J R) , 
with A = A + 2b~ 1 C 1 AQ 2 + C 2 b~ 1 h. Similarly, we have, 

N(Ah,Pv 2 (h)) < ■Af(^,b,Rh 1 ^ 2 ,R) . 

Thus we get from (|2.18p . 

N(Ah,P U2{h) ) < 2Af(X,b,Rh^ 2 ,R) . (2.19) 

Select ho sufficiently small and Aq sufficiently large such that 

2b- 1 C 1 A - 2 + C 2 b- l h < A - Ai . 

In this way, we get for all A £ (0, Ai) and h G (0, ho], A € (0, Ao]- Consequently, it follows from 
Lemma 12.41 that M( A , b, Rh 1 ^ 2 , R ) is bounded by a constant times (R 2 + 1), and thereby get 
the upper bound in (|2.17p . □ 

3. Asymptotic number of low-lying eigenvalues near corners 

3.1. The operator in an infinite sector. Let a G (0, 2ir) and 

= { (r cos 9, r sin <9)£M 2 :r>0, 0<<9<a}. (3.1) 

Define the magnetic potential, 

y( Xl ,x 2 )eR 2 , A (x 1 ,x 2 ) = (-x 2 ,0), (3.2) 

whose curl is constant and equal to 1. 

Consider the self-adjoint operator PQ a = —(V — iAo) 2 in L 2 (Q a ) whose domain is: 

D(Pn a ) = {u€H\n a ) : P Ua u G L 2 (n a ) , v ■ (V - iA )u = on 8fl a \ {0} } , 

where v is the unit outward normal vector of dVl a . 
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It is proved in [2] that the bottom of the essential spectrum of Pn a is a universal constant 
00 G (5) !)• Furthermore, when a G (0, the operator Pn a has discrete spectrum below 0o- 

Many questions connected with the spectrum of the operator Pn a are left open. Among these 
questions are the following: 

• Is the spectrum of Pn a below Oo finite ? 

• When a = ir, the operator Pq^ is the half-plane and its spectrum is purely essential and 
consists of the interval [Oq,oo). Is [0o,oo) the essential spectrum of Pn a for any a ? 

• It is conjectured in [5] that the ground state energy (bottom of the spectrum) of Pfi a 
is constant and equal to ©o when a G [ir,2ir). This conjecture has not been proved or 
disproved yet (see [2]). 

3.2. Assumptions on the domain. We describe precisely the regularity properties of the 
domain Q. The assumptions will be the same as those made in [21 E] . 

In this and the subsequent sections, Q is an open and connected set in M? whose boundary is 
compact and consists of a curvilinear polygon of class C 3 . By saying that the boundary T of 17 
is a curvilinear polygon (of class C 3 ) we mean the following (see |Grl p. 34-42]). For every x G F, 
there exists a neighborhood V of x in M 2 and a mapping ip from V to M 2 such that 

(1) ip is injective, 

(2) ip together with ip -1 (defined on ^(V)) belongs to the class C 3 , 

(3) llnV is either {y G : ip 2 (y) < 0), {j £ !! : ipi(y) < and ip 2 (y) < 0}, oi- 
ly G : ipiiy) < or ip 2 {y) < 0}, where ipj denotes the components of ip. 

The boundary T of $7 is a piecewise smooth curve. We work under the assumption that the 
boundary T consists of a finite number of smooth curves for k = 1, . . . , m. The family (T^) is 
the minimal family of curves making up the boundary T. The curve r^ + i follows according 
to a positive orientation, on each connected component of T. Let denotes the vertex which is 
the end point of r^. In a neighborhood of d£l, define a vector field v^. For each k, the vector 
is the unit normal a.e. on Y^. 

Let E be the set of vertices the domain $7 has. Suppose that S / and consists exactly of 
N vertices. This assumption certifies that the domain f2 does not have a smooth boundary but 
only a piecewise smooth boundary. When £ = 0, the boundary of O, is smooth and the operator 
in L 2 (0) is the one studied in [il] [10]. 

At each vertex G S, let a Sk denotes the angle between and measured towards the 

interior. 

3.3. Main result. For each angle a, recall the sectorial domain f2 Q introduced in (|3.ip . Let 
Pb,fi a = — (V — ibAo) 2 be the operator in L 2 {Q, a ) introduced in Section 13.11 (with Neumann 
boundary condition on the smooth part of the boundary of the sector Q a ). The number 

n(a,A;b) =tr(l ( _ 00 , A) (Pn a )) = £ dim (Ker (P Qa - Aid)) (3.3) 

A<A 

is finite for all a and A < Go- 
lf b = 1, we write, 

n(a,A) = n(a,A;b = 1) . (3.4) 
By a scaling argument, it is easy to see that, 

n(a, A; b) = n(a, A/b) . 

Let the domain Q be as described in Section \3. 21 Define the magnetic Schrodinger operator, 

Ph,n = -(hV-iA) 2 , in L 2 (n) , (3.5) 
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where A £ C 2 (0;IR 2 ) is the magnetic potential, h > is the semi- classical parameter and 
B = curl A is the magnetic field. The domain of the operator Ph,n is, 

D(P h>a ) = {ue L 2 (0) : (W - iA) k 6L 2 (fi), k = 1, 2, 

Vj ■ (/iV — iA)u = on T/t , k = 1, • • • , m} . 

Define two constants, 

6=infB(x), &' = inf B(x) . (3.6) 

The main result of this section is: 

Theorem 3.1. Let the constant Oo be as defined in (|2,5p . A € ( — oo, min(0o&', &)) and N(Xh) 
the number of eigenvalues of the operator Ph.n below Xh counting multiplicities. Suppose that the 
magnetic field B(x) is selected such that b > 0. There exists a positive number ho such that, for 
all h 6 (0, ho), the following equality holds true, 

AT(A/,) = gn(a Sfc ,^. (3.7) 

The number ho depends only on the angles a Sk and the domain O. 

Theorem 13.11 gives the exact number of low-lying eigenvalues of PhU corresponding to corners 
in the domain O. When A > Go, we will see in Section H] that corners no more affect the leading 
order term of N(\h). 



3.4. Proof of Theorem 13.11 

Upper bound. Let p be a positive constant satisfying < p < 1. We start by choosing a partition 
of unity Xkh introduced in [2| Prop. 11.2] satisfying: 

\Xk,h\ 2 = 1,1] iVxfc^l 2 < Ch~ 2p in R 2 , suppxfe,/, C B( Zj ,c k h p ) , 
k k 

with the choice of indices such that 

2^ = Sfc and = 1 for all A; = 1, 2, • • • , iV , 

If k {1, • • • , m.} and z& 30, then B{zk, ckh p ) n 30 = and c k = \ min(tan a Sfc , 1) , 
If k {1, ■ ■ ■ , to} and z k £ <90, then B(z k , c k 

B(z k ,c k h p ) n (90 is connected and c k = \ min(| tan a s J, 1) . 

Recall that £ is the set of vertices of the domain 0. If u is a function in the form domain of 
P h o, define, 



g(u) = / \(h\7 - iA)u\ 2 dx. 
Jn 



The following decomposition formula holds true for every function u in the form domain of the 
operator P ht n, 

III 



lib 

q( u ) = ^2<i(xk,h u ) + Y i(xk,hu)+ <i(xk,hu) - h 2 ^ / l v x 

k=l k>m j>m k 



k h\ 2 \u\ 2 dx. 



3'- 



Using the upper bound on |Vxfc,h| 2 , we get the lower bound, 

m ,. 

q(u) >^2q(Xk,hu) + ^2 <l{Xk,hu)+ ^ q{Xk,hu) - Ch 2 ~ 2p \u\ 2 dx. (3.8) 



k=l k>m k>rn 
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Let k £ {1, 2, • • • , m}. It is proved in [21 p. 252] that by performing a change of variable y = ip k {x) 
and a gauge transformation (defined by a function (ft k ), the following lower bound holds true, 

q(Xk, h u) > I ((1 - CbP - Ch 29 )\(hV - iB k A )v k \ 2 - C7/A" 2 V| 2 ) dy . (3.9) 
Here G (0, 1) is any constant, = B(^), Aq is the magnetic potential in (|3.2p and 

The optimal choice of p and # is /? = 3/8 and = 1/8. This produces the lower bound: 

q(Xk, h u)> f ((1 - Ch^KhV - iB k A )v k \ 2 - Ch 5 / 4 \v k \ 2 ) dy . (3.10) 

Similarly, when k {1, 2, • • • , m} and z k G (90, by applying a change of variable and a gauge 
transformation, we obtain the lower bound (see |llj). 

q(Xk, h u)> [ ((1 - Ch l ^)\(hV - iB k A )v k \ 2 - Ch 5 / A \v k \ 2 ) dy , (3.11) 
Ju h v J 

where Uh = {y = (yi, J/2) G 1R 2 : yi € (0, /i 3 / 8 ) and ?/2 £ (0, 00) }. By inserting the lower bounds 
in (pTLOj) and (f3~TT]) into QMS) we obtain, 

<z(«)>E / ((l-C^JKW-iSfcAoKp-Cfc 5 / 4 !^! 2 )^ (3.12) 
fc=i 

+ E / ((i-c^KW-iBfcAoKp-ch 5 / 4 !^! 2 )^ 

fe>m ^ 
z k £dQ 

+ E - c^ 5/4 / 



I -zxl 2 dx . 



k>m 



By the variational min-max principle, we deduce the following upper bound on the number 
N(\h), 



A + Ch 1 / 4 



N(Xh) < Vnk,- ^—j- , 

ti K (1 - C^/4) B ( Sfc ) 



+ 



E N(p Bk , Uh ,^±§^) + E ^(ift.Afc + C?^ 



k>m k>m 
z k £dfl z k gdfl 

Here PB k ,u h = —(hV — iBjAo) 2 is the operator in L 2 (Uh) with Neumann boundary condition 
at 2/2 = and Dirichlet condition elsewhere, and Ppn = —(h\7 — iA) 2 is the operator in L 2 (Q) 
with Dirichlet boundary condition. The spectrum of PB k ,u h starts at @oB k h and that of -P/f^ 
starts above bh (cf. |11|). Thus, if A < min(@o&',&) and h is selected sufficiently small, we get, 

N(p Bk , Uh ,^±§^)=N( Ph %Xh + Ch^)=0. 

Since the spectrum of the operator Pn a below Qq consists of isolated eigenvalues, we get for h 
sufficiently small that, 

X + Ch 1 '^ 
Sfc ' 1-ChWJ ~ "V" Sfc ' B(s fc ). 
This finishes the proof of the upper bound. □ 
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Lower bound. The proof of the lower bound is similar to the one in and uses a bracketing 
technique. Let P^q be the self adjoint realization of the restriction of the operator 011 
functions vanishing outside the set, 

N 

k=l 

By the variational min-max principle, the eigenvalues of Ph,n are larger than those of Py l: n- Thus, 

N(Xh) > N(Xh, P h ,n) ■ 

We will show next that there exist constants C > and ho > such that, for all h 6 (0, ho], we 
have, 

N / a - r/? 1 /4 \ 

N(Xh, P h n ) > V n ( a Sk , , \ , A . (3.13) 

If h is made sufficiently small, then 

A-C/1 1 / 4 \ / A 



""(l + ChV*)B(s k )J V Sfc 'B(s fc )^' 

and we get the lower bound we wish to prove. 

Derivation of (|3.13p is easy. Let q be the quadratic form associated to Ph,(i- 

Select an arbitrary function u in the form domain of Ph,n- A matching asymptotic upper 

bound to (pTLOj) is proved in gj p. 252]: 

[ |(V -iA)u\ 2 dx < [ ({l-Ch x/i )\(hV -iB k Ko) Vj \ 2 -Ch bli \v k f) dy , 

JB(s k ,h3/8) Jn ask v ' 

where = B(s&) and v k is obtained from u by the change of variables y = ijj k (x) times a gauge 
transformation. Summing over k, we get, 

q(u)> f {{l-Ch l ^)\{hV -iB k Ao)v k \ 2 -Ch^ 4 \v k \ 2 )dy. 

By the variational min-max principle, we deduce that the eigenvalues of the operator P^ q are 
larger than those of the direct sum of the operators Pu as ,B fc , thereby proving (j3. 13|) . □ 

4. Energy and number of eigenvalues in piecewise smooth domains 

4.1. Main results. We will state other results concerning the number and sum of eigenvalues 
of the operator P/^n introduced in (|3.5p . The assumptions on the domain Q is as described in 
Section 13.21 The notation and assumption on the magnetic field B = curl A is as in (|3.6p and 
Section 13.41 

Let (ej) be the increasing sequence of eigenvalues of the operator P^ q in the interval (—00, bh], 
counting multiplicities. If A E (—00,6], define, 



N(\h) = tr(^l ( _ 00iAh) (P^ Q )J = Card{e, : ej < Xh} , (4.1) 
E{\h) = tr(l ( _ 00i Ah)(i^,n " AMd)) = Y, ( e i " Xh ) ■ ( 4 - 2 ) 



The main result we prove in this section is Theorem 14.11 Its statement requires the notation 
[x]+ = max(x,0) for any real number x, the eigenvalue introduced in (|2.4|) . and the arc- 

length measure ds(x) along the boundary of Q. 
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Theorem 4.1. For any real number A < b and as h — > 0, it holds true that, 

IdQ J-oo 



B(x) 



dids{x) + o{h 1 ' 2 ) . (4.3) 



ul -2 P fOO 

E(Xh) = - h — / B(x) 3 / 2 
Furthermore, if A < b, then, 

N ( Xh ) = TTT/2 / B ( X ) V2 ^ ds ( x ) + i h ~ 1 ' 2 ) ■ ( 44 ) 

While proving Theorem 14.11 we only give the new ingredients and constructions required to 
adapt the proof given in |10| . We will refer to |10) for the detailed calculations. 



4.2. Proof of Theorem 14. li Notice that (|4.4p is a consequence of f|4. 3|) . Actually, the term on 
the right side of (|4.4p is the derivative with respect to Afo of that on the right side of ()4.3[) . On 
the other hand, using the variational principle of Lemma |2.2| the number N(Xh) can be seen as 
the derivative of the energy E(Xh) with respect to Xh. 

Therefore, to prove Theorem 14.11 it is sufficient to establish the asymptotic formula in (|4.3p . 

Using a magnetic Lieb-Thirring inequality and a decomposition of the operator by a partition 
of unity, it is possible to prove that the energy E(Xh) is finite for all A < b. Details regarding 
this proof are given in |10| Section 5.1, pp. 242-243]. 

We start by proving the asymptotic lower bound: 



Ll/2 f rCO 

E{Xh) > —-- — / / B(x) 3 / 2 



MO 



dids{x) + o{h 1/2 ). (4.5) 



B(x) 

Let H = Ph,n — Xhld and notice that the energy E(Xh) can be expressed in the more useful form, 

E(Xh)=tr(Hl { _ oo>0) (H)) . 

Let {/i, /21 • • • j In} be any L 2 orthonormal set in D(H). We will give a uniform lower bound to 

N 

Y. 0- ///; • 

Using Lemma |2,2| this will imply a lower bound to tr (HI^^ q^H)) . 
Consider a partition of unity of R, 

V ; i+V'2 = 1 i supp Vi C] - oo,l[, supp </> 2 C [-,oo[, (4.6) 

such that, 

|V^i| 2 + |VV 2 | 2 <C7, (4.7) 

and C > is a universal constant. 
Let 

r(h) = ^min(| tan a|,l)/i 3/8 , (4.8) 

where a = min{a Sfc : k = 1, 2, ■ • • , m}. 

Using the partition of unity in (|4.6p . we put 



t(x) = di S t(x,dn), ^ h ( x ) = ^ 1 (J^±\ <M*) = ^ (^) , Vsen. (4.9) 

We introduce the potential 

V h = h 2 (\V^ h \ 2 + \V^ h \ 2 ) . (4.10) 
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It is possible to prove that (see |10|. Proof of (5.10), pp. 243]): 

N N h / h \ 

E<£ , Hfj) > , ^ h {H - V h )^, h fj) - C— \l + -^j . (4.11) 

With our choice of r(h), the remainder terms in (|4,lip are of the order of h 5 ^ 8 = o(fo 1//2 ). Thus, 
we have, 

N N 

E</j , Hfj) > Yfrj > M H - V h)^i,hfj) + oih 1 / 2 ) . (4.12) 

3=1 3=1 

Next, we estimate from below the term on the right side of (|4.12p . Recall the vertices Sfc, 
k = 1, 2, • • • , m, of the domain Q. Recall also that the boundary of the domain £1 consists of 
smooth curves (1^). For each k, define, 

T k (h) = {xeU : dist(x,T fc ) < r(h) and dist(x, s k ) > ^ 3/8 } , (4.13) 

where r(h) is defined in (14.811 . Consider a partition of unity of ' 



p2 . 



TtL TtL TTL 

£^+5>2 = i, E(i v ^i 2 + iv^i 2 )<^, 

k=l k=l k=l 

where suppg k C B(s k ,h 3 ^ 8 ), supp/i/% C T k (h) for all k, and C is a universal constant. 
Then, we have the decomposition formula, 

N N m 

^2(fj , i>i,h{H - Vh)ipi,hfj) = E E^' ' 9ki>i,h{H - V lth )g k iJ}i,hfj) 
3=1 3=1 k=i 

N m 

+ E E^i ' h kM H - Vi,h)hk^i >h fj) > ( 4 - 14 ) 

3=1 fc=l 

where 

m m 

V hh = V h + h 2 J2 |V 5 fc| 2 + h 2 E |V/i fc | 2 , 

k=l k=l 

and 

\V hh \ < Ch 5 / 4 . 

We will show that the first term on the right side of (|4,14p is of order o{h 1 / 2 ). Recall the definition 
of the energy S COTn (a, b, R, A) in (|2.1ip . We use (|3.1ip with u = ipi,hfj and Xk,h = 9k to bound 
the term {fj , g k ipi^(H — Vi^hjOk^i^hfj) from below by a quadratic form defined over the sector 
f2 Qs . Then we use the variational principle of Lemma 12.21 and the lower bound of to get, 

N 

£</j > 9kM H ~ Vi,h)9k^i,hfj) > -(1 - C/i 1 / 4 ) <5 C orn(as fe , B fe , /i -1 / 8 , A) . 

3=1 

Here A = Aft" 1 - 1 + C/i 1 / 4 < Ch 1 / 4 . Clearly, A can be made smaller than an arbitrary positive 
number. Thus, we can apply Lemma 12.51 and use (|2.14p . In this way, we get the following lower 
bound that holds uniformly with respect to k, 

N 

E)(/i • 9kM H - Vi,h)9kA,hf3) > -Ch 3 ? 4 = oih 1 ' 2 ). (4.15) 

3=1 
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The second term on the right side of (|4.14p is bounded as in |10^ Proof of (5.26); pp. 244-248]. 
The following lower bound holds uniformly in k, N and the orthonormal family {fj}: 

N 



^2ifj > h k ip lth (H - V ljh )h k tp 1:h fj) > 



ft 1 / 2 
' 2vr 



B{xf' 2 



A 



B{x) 



d^dx - Ch 5/8 . 

(4.16) 

Substitution of f|4.15j) and (|4.16|) into ()4.14p establishes the lower bound (|4.5p . 

Upper bound. The upper bound will be obtained by constructing a specific density matrix 7 
and computing the energy of truP/^n — AId)7j. The calculations follow closely those in )10| 



Section 5.3, pp. 248]. 

Recall the sets Tk(h) introduced in (|4.13|) . Let r(h) be as defined in (|4.8p . For each k, we 
cover Tk(h) by disjoint squares, 

Kj,k = {x eM 2 : \x- zj,k\ <t}, l<j<N k , 

where the points (zj k ) are on T k and equally spaced. In each Kj ik , it is possible to apply a 
transformation <§>j k {x) = (s,t) that straightens dfl (t = defines the boundary of f2 and s 
measures the arc- length along d£l). 

Recall the eigenvalues (/^>(£) of the harmonic oscillator in (|2.3p . We introduce an orthonorml 
basis (u p ) of L 2 (R + ) : 



(4.17) 



Let x ^ C°°(M 2 ) be positive, smooth, have values in [0,1], supported in B(0,2) and equal to 
1 on -B(0, 1). Define Xj,k(x) to be 

' x Zj,k N 



4j0;£) =0, 



r(h) 



We denote by Bj ik = B(z J - ) fc). For each (j, k), there exists a function 0^ such that 
A - V^fc = f - x 2 , 0) + O (\x - zj,k\) in Kj,k ■ 

Define 

±5, / x 



/j,p,fc((s,^);/i, ,0 := 



We get — by the coordinate transformation — the following function in Q, 

fj, P ,k( x 'i h ' = fj, P ,k(®j,k( x )'i h ' Xj,k{ x )^i,h{x) ■ 
Next, let if > and define the function, 

m N k 
fc=l j=l 



(4.18) 



(4.19) 



where M p (h,£, j,k ; K) = M (h, k)lr p= x\ and M (h,£,j,k) is the characteristic function 
of the set: 



{&j,k) : 



A 



B 



mi(6>o, lei <x>. 
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Let Hp (h,£) be the operator with integral kernel in L (Q) whose integral kernel is defined by, 

U h p nd (h, 0(x, x') = f p (x ; f p (x';0 . (4.20) 

Define the density matrix, 

1 oo 

^^W^r^). (4.2!) 

Clearly < 7 as an operator on L 2 (il). Furthermore, there exists a constant C > such that, 

7 <(l + Cr(/0), (4.22) 

in the quadratic form sense. Details concerning the derivation of (|4.22p are give in |10^ Proof of 
(5.34)]. Moreover, (see [TUJ (5.37)]), 

tv[(P h!U -\hh}<-^ [ K [ B(xf/ 2 [-±--^)} + ds(x)dt + CKh 3 /*. (4.23) 
2vr J-KJdn L B(a;) J + 

Let us mention that, while estimating the term on the left side of (|4,23p . the discrete sum over 
j and k becomes a Riemann sum as h — > 0, thereby resulting in an integral on the right side of 
(|4.23p . Since K can be chosen arbitrarily large, (|4.23p implies the asymptotic upper bound, 

l1/2 roo r \ 

tr[(P ht a-Xhh]<-—- / / B(xf 2 [—-^)] + ds(x)dti + o(h 1 / 2 ). 
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